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ABSTRACT 

It is pointed out that insight into the threshold region of 
electron -atom ionization can be gained by examination of the nature 
of the doubly excited states of the compound ion. A study of these 
states for H~ has been initiated with two types of variational wave 
functions. One, Y^, has the two electrons at roughly equal distances 
from the nucleus; the other, has one electron at a very much 
farther distance such that it sees the dipole potential caused by the 
inner electron and the nucleus. Both functions are constructed to be 
eigenfunctions of the operator Q^, which projects out all states of 
the target of principal quantum number less than N, and renders the 
energy subject to a minimum principal. If the number of states 
for which Yy yields a lower energy then y d is proportional to N^, then 
an extrapolation argument shows that the threshold yield curve 
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will be effectively proportional to Calculations 

have been done to N - 5 for Y n an d N — 9 for Y • Only the 

D W 

lowest of the Yy states is lower than the corresponding Y^ state. 

The results suggest y in the range 0 < -y s - l/2. Some comments 
on Wannier 1 s theory of ionization are made. 

I. INTRQDU CTION 

The theory of low energy electron impact ionization of atoms 
by electrons is fraught with difficulties from begining to end. 

The difficulties are both mathematical and conceptual in nature. 

The mathematical difficulties derive from the long range nature of 
the Coulomb potential combined with the intrinsically three -body 
nature of the wave function in the final state. In almost all 
cases, however, these problems are related to conceptual questions 
of immediate physical significance. If the two electrons come away 
from the nucleus or residual ion (considered an infinitely heavy 
point charge and always referred to as the nucleus) with approximately 
equal and opposite velocities, then it is a reasonable argument that 
each electron sees the nucleus directly and that classical mechanics 
can be applied 1 . The point here, of course, is the virtual identity 
of Coulomb scattering in classical and quantum mechanics and the fact that 
the classical approximation becomes more exact as the energy gets lower. 

If on the other hand the electrons come off with quite different 
velocities, the validity of classical mechanics is a much more 
questionable item. Here the quantum mechanical argument is that the 
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inner electron may shield the outer electron from the nucleus 
thereby making the potential it sees shorter range than Coulombic. 

In that case its behavior may not be governed by classical 
mechanics and the classical approximation may get poorer as the 
energy is lowered 2 . 

It is the purpose of this paper to examine some of these 
questions from a consistently quantum mechanical point of view. 

We shall attempt to avoid questions concerned with the controversial? 
asymptotic form above threshold 3 * 4 by considering the process as 
a continuation of real or virtual processes below threshold. It 
also allows calculations to be done in a fairly unambiguous way. 

We believe that this is the most important aspect of our work, for 
if any question be raised concerning the variational forms of our 
wave function^ the way has been opened for other forms to be proposed 
and tested on the impartial balance of quantitative comparison. The 
one question that this approach can probably not answer is any subtle 
questions of analytic continuation from negative to positive energies. 
We shall be more detailed concerning what effects we believe this can 
have in Section III; for the present it is only relevant to note that 
the one advantage of the Coulomb force is that the continuum solution 
merge continuously with the discrete solutions. Thus if we correctly 
describe the major physical situations that can occur below threshold, 
then we can be reasonably sure that what we extrapolate them to be 
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must be substantially correct. 

What are the physical processes which extrapolate to 
ionization? They are of two types : (l) Inelastic processes in 

which the orbital electron is raised to a highly excited state N 
with N finally going into continuum. (2) Double excitations in 
which both particles can be considered simultaneously caught in an 
excited "bound" state, also characterized by the principal quantum 
number N, where again N finally goes into the continuum. Whereas 
the former are real both in the sense that they occur at energies 
above the energy necessary to excite the Nth level and the processes 
correspond to rigorous time independent solutions of the SchrfJdinger 
equation, the double excitations processes can be considered virtual 
in as much as they occur below threshold and they do not correspond 
to rigorous time independent solutions of the Schrbdinger equation 5 . 
Nevertheless they do occur, and they can have a profound effect on the 
scattering both below and above the thresholds in question. 

Our approach will be to use the double excitation processes to 
guide us in the choice of final state wave functions that we use, 
but the actual derivation of the threshold law will be carried out 
via (l). In Section II we deal with doubly excited state calculations 
and in Section III with the derivation of the threshold. One of the 
forms of the doubly excited wave function has been motivated by 
Wannier 1 s theory of ionization 1 according to which the threshold law 
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is dominated by processes in which the energies and the radial 
distance of the escaping electrons are not too different from each 
other. In Section IV we discuss his theory a little more and point 
out that his derivation cannot be justified in a completely classical 
theory of electron -atom (or ion) ionization. 

It turns out that the scattered particle must see at the very 
least an r~ s potential. This causes a change in the continuum 
normalization factor of the scattered wave which makes it more like 
a Coulomb wave. This normalization factor is derived in the appendix; 
but it is also indicated there that the normalization factor notwith- 
standing, the ionization in such a potential should have an energy 
power dependence less dominant than that of a pure Coulomb wave. 

II. DOUBLY EXCITED STATES OF H - 

We shall calculate doubly excited (i.e., autoionization) states 
of H~. There are two reasons for dealing with this negative ion: 
first it is the simplest negative ion, but more important the eigen- 
functions of the target atom are knowi exactly and therefore we can use 
the Q-operator technique 6 ’ 7 without approximation. Aleo following 
Wannier 1 we s hall deal with only total S-states in the belief that 
threshold law cannot be altered in form by higher angular momentum 


states . 
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The two wave functions which we use are?: 




R 4 (r ) 
vi 1 


'vNjf 



v :> N 


( 2 . 1 ) 


and 







(-i) x /2 Tl T p^cos e 12 ) 


( 2 . 2 ) 


(n) 

The physical meaning of y is easily understood, it describes 

the two electrons in doubly excited states at roughly equal radial 

distances from the nucleus. (We asume v = N). The functions Rj^(r) 

are r times the radial hydrogen wave function. The angular 

correlations indicated by the electron-electron repulsion are taken 

up by the linear variational parameters C „ of which N - 1 are 

vlU 

effectively free and one is determined by normalization. Physically 
it is clear that the calculation will make them such as to concentrate 
the electrons on opposite sides of the nucleus. It is also clear 
that for neither particle does this function contain states of 
hydrogen with principal quantum less than N. I.e, for i _ 1, 2 
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J % (r i> V n i> v <N)d3r i = 0 


(2.3) 


n < N 

\ <z £, <; N - 1 


Thus the function is an eigenfunction of the operator 9 and 
variational calculations will give eigenvalues which if they lie 
below the Nth state of the hydrogen atom will correspond to 
resonances in the elastic and inelastic channels. 

Although the physics of is also readily understood, the 

mathematics needs some explanation. The following is a precis of 
Mittleman' s 1 ° generalization to arbitrary N of the analysis which 
Temkin and Walker 11 have given for the N _ 2 state. Let us start 
with the following ansatz for the closed channel wave function : 


f o (N)= A W r2) + 1 ' P / cos O . 

I = O r n r 2 


( 2 . 4 ) 


The functions u^(r) are to begin with undetermined functions. [The 
factor (-1) V 2 j& + l) is the essential part of the Clebsch-Gordan 
coefficient by which Y^fa,) andY^ _ m (n 2 ) 00uple tQ form p , )y 

AJ 

If one varies the in the expression for expectation of the energy, 
then one arrives at coupled differential equations whose longest 
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range terms are of the order r~ 2 and to that order the equations 
may be written 


(N) I (N) 


dr 2 r 2 


U 


(2.5) 


(n) (nI 

U is a column vector of the and J3 is a tridiagonal 
symmetric matrix whose elements are given by: 


-(H) Jt(i + 1) r ~,.i A 2 - l' z ' . 

V 3 ~~r^ hi + 5N * 7 . x h + 1. i 


5NA 


/N 2 - ~E 

bf - 1 


6 A, A' + 1 


( 2 . 6 ) 


Introducing the transformation 

V (N) _ b (N) u (N) (2.7) 

such that is diagonal, we find that the components of 

(n) 

the column vector V satisfy the equation 


r_sL - !si 

L dr 2 r 2 


N, 


+ = J » 1) j 


( 2 . 8 ) 
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(n) (n) 

where b„ T ,are the eigenvalues of 45 . Let B ' be the associated 

WJ jW 

eigenvectors, then the function Y^ of Eq. (2.4) takes on the 

form of Eq. (2.2) and acquires an additional "quantum" number j which orders 

(n) 

the sequence of eigenvalues of £5 and their corresponding eigenvectors. 

Of the N eigenvalues b„. a certain number, j _ 1, 2, — , J , 

it J J) 

will be negative and for those functions the equations will asymptotically 
contain an attractive r 2 potential. For each N and j the equations 
will contain an infinite number of negative eigenvalues which to 
an excellent approximations are related by 1 © 


I U + Di 

e Nj 


- ^>1; s 


1 , 2 , • • • , 


(2.9) 


where 




[|b BJ | - iM 1/2 


( 2 . 10 ) 


The solutions of those equations are Hankel functions 11 > 12 

H. (ilewJ^r) which asymptotically approach exp(-| e „, | 1//2 r) . 
la Nj NJ Nj 

At short distances equation (2.8) becomes altered and essentially 
non-local in character. In fact the attractive r~ 2 potential must 
become less singular, for the r" 2 solutions are not regularly behaved 
at the origin. 
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For variational purposes we have taken 

<1 

v (r) = e” ar V C r m . ( 2 .ll) 

NJ L m 

m = 1 

In order that he an eigenfunction of it is necessary that 

00 

I V r) R n X (r) = 0 < 

o 

Note that the Qjj operator is symmetric 7 • Although the form of 

as it stands appears to contain no bound states of lower N 
for the target particle (r 2 ), it might contain lower states in the 
exchanged coordinate (r^ ) . In order that the calculation be 
subject to a minimum principle the possibilities of ordinary 
exchange inelastic scattering from a state lower than N must also 
be excluded. Eqs . (2.12) @iarantees this to be the case. They 
are a set of q - 1 linearly independent equations where 


n < N 
\ < N - 1 


( 2 . 12 ) 


q. = 1 + N(N - l)/2 . 


(2.15) 


If the v ( r ) contains exactly q terms, then together with normal- 
ization all the coefficients are unique functions of a. Thus 
for variational purposes the function contains only one 
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(n) 

variational parameter, a, as opposed to ' which, as was stated, 
contains N - 1 parameters. It is very important to realize, because 
the W and D calculations are based on the same operator, that the 
shifts Aq are the same and therefore can be omitted in the comparison 
of the respective energies that will be made (Table III) . 

(N) 

The matrices B were inverted for all H = 2 to H = 100. 

The number of negative eigenvalues is clearly linear with N as 
is evident from Figure 1 in which is plotted as a function 
of U. In fact it is quite certain that 

= \ ^ " (lower order in N), (2.14) 


where the term in brackets may very well be logarithmic. The 
eigenvalues themselves appear to go up quadratically with N. 
Selected values are given in Table I. The differences are seen 
to be proportional to N and independent of j for the lower values. 
In fact the b„ can be fit to a reasonable approximation by a 

J 

formula of the form: 


b ^ 3 N 2 + B(j) • N + C(j) (2.15) 

Nj 

(n) 

Some feel for the eigenvectors B^' can be gleaned from 
Table II. There it can be seen that aside from a normalization 
constant N 1 / 2 the eigenvectors are fairly independent of j for 
small Z and that they get extremely small fur large » as long as 
j is small. For large j they lend to oscillate and they are all 
of the same order of magnitude. 
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The significance of these properties will emerge in next 
sections in which we utilize these properties to derive an effective 
threshold law. We shall also discuss there the significance of the 
results of the variational calculations. The results themselves 
are summarized in Table III. 

III. IMPLICATIONS FOR THRESHOLD IONIZATION 

It may appear that the variational wave functions we have 
used, in particular Yy> are overly restrictive and that specifically 
if we had used shielded Coulomb radial wave functions we would have 
found more energies of the two Coulomb type (Ey) lower than their 
dipole (Ep) counterparts. (This number is indicated by Jy in 
Table III.) We do not believe this to be the case for the following 
reasons. The number of variational parameters in Yy goes up as N, 
therefore one has more freedom to simulate the effects of shielding 
should this have been required. To test this point even further 
we have varied v in Yy away from N. The effect of this change can be 
described as making the mean radial distances of the electrons be 
slightly different from each other, and thus to give partial shielding 
more room in which to operate. The results are summarized in Table IV. 
It can be seen that in only one case does it reduce an eigenvalue 
(v - 10, N = 9 ) and there only the lowest one (which was lower 
than Ep anyhow) . Its effect on all higher eigenvalues is to raise 
them, and in fact in only one case does the second eigenvalue remain 


bound 
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Furthermore it must be added that the function y (for v __N) 

w 

is explicitly symmetric whereas y^ is not. Although we do not 
expect this to be a significant factor for large N, Table III 
shows that it can have a significant lowering effect for small N. 
Finally we reiterate that y^ only depends on one variational 
parameter. Thus to the extent that this calculation is biased, 
it predominantly favors the two Coulomb functions y^. This is 
intentionally done to offset any criticism that we were intuitively 
drawn to the shielding approximation 2 . 

Thus it is the calculations themselves which convincingly 
demonstrate that the y^ functions dominate. In other words in 
double excitation the electrons tend to be at greatly dissimilar 
distances from the nucleus. For example the second autoionization 
state of the N - 5 calculation has a mean radius of r 2 at approxi- 
mately 37 a Q whereas r 1 "1 136 a Q . This in turn can be described 
as the direct effect of shielding of the outer electron from the 
nucleus by the inner electron 2 . 

The threshold law for ionization, however, will be determined 
by those few states in -which the electrons emerge at comparable 
distances from the nucleus. The present calculation clearly shows 
that such equal energy events can occur . In order to extrapolate a 
threshold law from the present results we shall proceed as follows 13 . 
To every two-Coulomb autoionization state below the Nth threshold 
we shall associate cm. xxicXsl stic scattering wave function above the 
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Nth threshold in which the inelastically scattered particle will 
also be described by a Coulomb wave. This wave function may be 
written 


*y. > J o 

L N&j 

i = o 


r 2 


(cos e 12 ) 


(3.1) 


The index j here labels the state of the calculation whose 
energy is lower than the corresponding j state of the y d calculation. 
In general then 


J = 1, 2, ..., J w (3-2) 

Although for the results presented in Table III J^. £ 1, it must 
be assumed that as N gets larger, will also. We shall assume 
that asymptotically for large N, can be represented as 

11m oc N Y , 0< v sl 

N-co W (3.3a) 

The function however, is an inelastic scattering wave 

function, and is wave number of the inelastically scattered 
wave. The cross section for excitation of the Nth state is 


given by 
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j = 1 


< I v I ** > I 2 


where is the initial state and V is the interaction: 


§. = t cp (r 2 ) ] 


sin kr x R 12 (r 2 ) 


<Po K±2 '* TEr” - 

S wave i 


2 2 
V = ~i 


(Rydberg units are used throughout). 

We define the yield to a group of states in the vicinity of the 
Nth state as Q (not to be confused with the projection operator 


N + AN 

Q = Y a N 

N - AN 


Now if the energy of the hydrogen atom is labelled w(N), then 


w(N) = - — 
N 2 


(3.8a) 



- 16 - 


the total energy being given by 


E = l ^ 2 + w 


( 5 . 8 b) 


When N is large, we can replace the sum in (3-7) by an integral: 


Q 


ajjdll oc J ctjjN 3 dw 


(3.9) 


The matrix element in (3*4) can be written explicitly 
00 

<%jl v I »i > “I J J dr > dr2 p( v. )( w V r > 1 (jiTT> x 

i = o 


V r2> 


1 + 1 


®is ^ r z) 


( 3 . 10 a) 


14 - 3/2 

The normalization factor of is proportional to N . 


Since there are N terms in y. T . the coefficients 


like the B 


(N) 


Xjj uie coeumieuiia - L - LJ% - C otlc “ $ * 

- l /2 

Table II, must be proportional to N ' . The normalization of 

14 

F(kjj r ! ) at the origin, corresponding to unit current at infinity , 

- 1 /2 13 
is kjj ; the sum of integrals over 1 converges very rapidly 

and is quite independent of the upper limit. We therefore have 




v 1 *, > « — r 


1 N2 


( 3 . 10 b) 
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We shall assume that the dependence on j for N large is 
secondary, as with the coefficients b^ of Eq. ( 2 . 15 } so that 
the sum over j in (3*4) contributes 15 a factor J^. Eq. (3-9) 
then leads to 


Q OC 



f J w N 3 dw 


(3.H) 


Using ( 3 . 3 ) for J^, inverting (3-8a) in the form N a | w | 1 / 2 

and proceeding into the continuum wherein w is positive and 
0 £ w * E defines the range of integration, we find that 



(3-12) 


Thus finally 


Q oc E 2 (3-13a) 

In order to say something about the value of y, we have 
plotted in Fig. 2 the results we have obtained as a function of N. 

The solid straight line is the same as that in Fig. 1 in the 
restricted range of N. It represents and therefore is an average 
of the squares; although it looks somewhat arbitrarily drawn here. 
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Fig. 1 shows that when one goes to larger N there is essentially 
no ambiguity in it. Similarly the curved lines are intended to 
represent an average through the open circles (J.^. vs. N). There 
are not enough open circles to allow such a curve unambiguously 
to be drawn, however as with we expect that going to larger N will 
allow this curve to be essentially uniquely continued. The limited 
results do seem to fit better with the smaller fractional value of 
y. Further discussion is reserved for the next section. 

IV. DISCUSSION 

Let us examine some of the assumptions that have gone into the 
derivation of Eqs. (3.15a). Aside from the analytic continuation 
into the positive energy domain, the biggest assumption concerns the 
summation over j in (3 •*<•)• We have indicated in footnote 15 that 
the only error this could reasonably cause is an increase in the 

exponent in (3 . 13 a) • To that extent the exponent -X 

V 2 

be a lower bound on the exponent, which would be quite satisfactory 
for our purposes. We believe, however, that it is more accurate 
than that. The process of analytic continuation which is used in going 
from (3.1l) to (3*12) together with the restricted analytic form used 
to represent may have lost more subtle energy dependent factors 

such as log E or oscillating factors. From most practical points 
of view logarithmic terms are not important, since they are completely 
dwarfed by the power dependent factors; if the factors are oaoillating, 
then we would expect (3-13) to describe the envelope of the curve. 
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The smallness of in the present calculation might he in interpreted to 
mean either that y is small or that the constant of proportionality is 
small in (3.3a). We believe that it is the former which is more suggested by 
our results, for otherwise one wouldn* t have expected J^.(N) to be 1 at N _ 3 
and then to remain there to at least N _ 7« Rather one would have expected 
only to become 1 at a larger value of N. It is for this reason that the 

agreement of with a fractional power dependence in Figure 2 appears to fit the 

limited results so natural lv. Note also from Table III that J TI cannot exceed 

W 

2 for N <: 9. 

Nevertheless it will clearly require much larger N in order for a precise 

value of y to be determined. Although we axe in the process of extending 

this calculation, we can not promise that results will be forthcoming soon. 

It will require considerably more numerical sophistication to avoid overflow 

and cancellation of significant figures. (The computer is an IBM 360-91 with 

approximately 15 significant figure accuracy. ) 

On the basis of Figure 2 we would estimate 0 < y ^ l/2. Indeed the most 

likely alternate possibility in our opinion would be a logarithmic increase 

of J Tr with N, 

W 

J w cc log N (3.3b) 

This would change the form of the threshold law to: 


3/2 

Q oc Yr' log E 


(3.15b) 


Although we think it is unlikely, our results are not extensive enough 
to rule out Wannier' s 1 threshold law. (This would correspond to y = 0.75* 
A linear theory 4 y = 1 seems distinctly improbable in our opinion) . This 

X 6 

theory, "which has recently been revived by Vinkalns and Gaiiitis 
is based on a rather brilliant analysis of the classical 


9 



- 20 - 


orbits (i.e. solutions of Newton's equations) which describe two 
electrons emerging from the vicinity of the nucleus and not being 
caught again. Basically Wannier finds for E zero or slightly 
greater than zero that the solutions are of two kinds: 


Ar 


C, r i/k - " /2 


(4.1a) 


Ar = C 2 


r 5A + n/2 

r 9 


(4.1b) 


where &r is the difference of the radial distances of the two 
electrons fran the nucleus (assumed fixed and of charge Z) and 
r is the mean radius : 




+ r 2 e 


(4.2a) 



(4.2b) 


and 


M- = 



100Z - 9 
4Z - 1 


(4.3) 


The first type of solution (4.1a) can exist even at E = 0, and it 
corresponds to particles appearing at infinity with equal (necessarily 
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zero for E = 0) speeds. Geometrically similar solutions continue 
to exist for E > 0, thus the threshold dependence (increase in 
the number of solutions) is determined by the increase with E 
of solutions of the second kind (4.1b). These correspond to events 
in which th6 two electrons come off with slightly different energies. 
Wannier 1 has shown that the contribution of these to Q is proportional 
to C 2 itself. To get the dependence of C 2 on'E Wannier appeals to 
a similarity principle whereby if 

r = r (t) , i = 1, 2 (4.4a) 

~uL 

is a solution of Newton 1 s equations for energy E, then 

r ' = h r. (b 5/ 2 t) (4.4b) 

~1 ij 

are (geometrically similar) solutions for energy E 1 - BE. The 
solutions (4.1b) can be written as an explicit function of time 
using the solution of joint motion 

r « t 2 ^, (4.5) 

which is valid when E « Zr \ This can then be consistently 
inserted in (4.1b) (which itself is valid only when a r « r.). 



- 22 - 


Letting C 2 inaX (E) be the maximum value of C 2 which leads to "double 
escape" at energy E, one finds that the corresponding solutions 
(4.1b) can be written: 

. *<»“>(.) ^ + ^ fc. le 

Applying (4.4b) then leads to the conclusion 

C 2 inax (E) oc E “ Z + 2 (4.6) 

In Eq. (4.6) the (qjiasi ergodic) assumption has been made that all 
initial conditions for particles entering the emergent zone are 
essentially equally probable. From the remark above Eq. (4.4a) this 
then translates itself into 

Q « E~ + ^/ 2 , (4.7) 

which is Wannier* s threshold law 1 . 

The key assumption in this theory, in our opinion, is expressed 
in Eq. (4.6). We wish to show first that this assumption cannot be 
Justified in a strictly classical theory of the whole ionization 
process. In that case the cross section emerges as a statistical 
average of events in which the orbital particle is initially bound. 
(Radiation damping is necessarily excluded). The question one asks 
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in ionization is what happens as the energy of the impinging particle 
increases, the characteristics of the bound particles remaining the 
same. In other words the variation with energy of -the initial condition 
with energy does not satisfy the similarity principle. Eqs. (4.4). But 
Newton' s equations in that case cover the whole collision process . In 
other words if the orbits corresponding to the solutions (4.1) be 
traced backward in time, it will be found that the overwhelming majority 
of them originate in trajectories in which the two electrons were 
originally approaching the nucleus from infinity. These are initial 
conditions that must be excluded even from the most general type of 
distribution used to describe the real initial conditions. Thus we 
conclude that from a completely classical, point of view the distribution 
of C 2 ^ maX ^(E) does not necessarily obey (4.6). It may be a very 
very sensitive function of E (near threshold) and/or the result may 
not depend on C 2 (e) alone but on C x as well; it may also depend on 
the statistical distribution that one chooses to describe the initially 
bound orbits. It is, then, perhaps significant that classical monte-carlo 
calculations, do not reveal the Wannier threshold law. 

The same objection cannot a priori be raised against a quahtum 
mechanical collision. For in that case the concept of an individual 
orbit does not apply throughout the collision process. Nevertheless 
the above consideration does raise the likelihood that the probability 
with which particles emerge into the classical zone may also be a 
highly sensitive function of E. We believe that this in fact is the case. 
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since even if one cannot speak in terms of orbits, the Ham iltonian 
does remain the same throughout the collision. The assumption 
would be more justified in our opinion in the problem of the 
threshold production, say, of two negatively charged hadrons by 
nucleon collisions with nuclei. In that case the short range in- 
teractions only come into operation in the quantum mechanical zone, 
and they are so strong and complicated that they can legitimately 
be expected to make the final state completely oblivious to the 
initial state. 

In addition to this quasi-ergodic assumption there remains the 
question of the validity of the classical theory. This is a very 
difficult question which has not been definitively answered. We 
believe that in the ionization of atoms by electrons the theory does 
have some validity in the region (&r/r) « 1. However when the 
difference in the two radial distances gets large, then we believe that 
quantum effects (shielding) will have a profound role 2 . The present 
results tend to bear out this reservation. 

Finally we mention the experimental situation. Although an 

experiment can never prove a threshold law, the experiment of McGowan 
is 

and Clarke has convincingly shewn that there is some nonlinearity 
in the e - H ionization curve between threshold and 0.4 eV. For if 
there were not then the measured position of the first resonance in 
e - H elastic scattering 


would not coincide with essentially pre- 
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cision calculations, which in our opinion can not seriously be 
questioned. In the region 0.05 < E < 0.4 eV McGowan and Clarke 
find very good agreement with Wannier 1 s law, but it is 

perhaps (perhaps because the experiment is hardest there) significant 
that below 0.05 eV the yield curve does appear more nonlinear. Brion 
and Thomas 20 in e - He ionization also find a yield curve which 
appears to be more nonlinear than Wannier 1 s law. Since the region 
of nonlinearity is much greater there, this may provide a better 
experimental test of the theory. 

Appendix 

In this appendix we shall derive the normalization factor for 
a particle scattered in an attractive r 2 potential. We shall also 
estimate its effect on the ionization threshold. 

The inelastically scattered wave satisfies the equation 



(Al) 


It is important that this equation exclude regions near the origin; 

21 

this is evident from the general solution of (Al) 


f (r) = Ar l/2 J 
Nj 1Q?| 


(V) + Br ' 


1/2. 


N 


Nj 


i'i 


ia, 


(r). 


'NJ 


(A2) 
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which oscillates infinitely rapidly (and thus is unacceptable) 


as r -> 0. 




is given by 





(A3) 


Using the asymptotic form of the Bessel function for k^r » 


in (A2) gives 


m 

i - . V r) = A /3^ cos ^ ■ f ,+ B sln (v • ^ • £> 


We demand that f^(r) be normalized to unit current, so that 


llm fgj(r) = j- sin (y + 6 ) 

r ^ co N 


(A5) 


where C is independent for k^. Comparison of (A5) and (A4) shows 

-l/2 

that A and B are proportional to k^ ' . 

On the other hand A and B must also be related to the solution 
for smaller vhlues of r. To get this relation we must first rewrite 
(Ah) in a region -where k^r « q^j, but where the r 2 potential is still 
operative. The function f^(r) can there be approximated by the small 
argument expansion 22 of the Bessel functions. I.e., for kj^r « o^j 
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but r > r Q 


(A 6) 


f (r) ~ Ar 1 ^ 2 C ^ ) Br 1 / 2 r/'V\ * cos (iorr) “ *1 

NJ =Ar fa ~ia f sin (I^TV rll + X) ' fCl - 1«) J 


The radius at which the r" 2 potential is no longer operative is 
defined as r Q . Within this radius we assume that we can represent 
the solution for whatever (in reality very complicated and nonlocal) 
potential does exist as 


f H j(r) = T)f(r) 


(A7) 


The quantity T) is desired normalization factor. It can be determined 
by equating the logarithmic derivative of (A6) and (A7) ar r _ r Q . 

One finds to an excellent approximation 


1/2 


T1 = 


f(r o> - Vo + R oV ) 


(A8) 


where 



= r 

o 


r 


(A9) 
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In order to estimate the effect of this on the ionization 
it is necessary to have some idea of the size of r Q . A very 
reasonable estimate of it is that radius at which the outer electron 
is comparable to the mean radius of the inner electron r 2 . (it cannot 
be smaller for then it would no longer be the outer electron.) But 
for a hydrogenic atom in the Nth state 7 S a N 2 , thus it is clear that 
as N gets large r Q must also, hence 


lim 
N - 00 


11 


' f <ro)| ( k N r o) 


1/2 


(A10) 


Before proceeding let us note that for a finite N the factor 
in (A8) or (AlO) is the same as one would have for a pure Coulomb 
wave. It is this factor which is directly responsible for the finite 
threshold behavior in electron impact excitation of hydrogen. This 
simple fact which is the implicit basis of the original derivation 
of this result by Gailitis and Damburg 21 is somewhat obscure 
in their paper as a result of their very elegant and very general 
mathematical procedure. 

To calculate the effect on the ionization we must replace F(k^r 1 ) 
in (3*1) by T]f (r 1 ) in view of the fact that the dominant contribution 
to the matrix element comes from r 2 small and r x comparable to r 2 . 

In this case the sum (3-4) must be extended from j _ 1 to j - J-p. 


preceding 


page BLM« not BIMEO. 
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Figure Captions 


Figure 1. 


Figure 2. 


J D vs. N. Both and. N are defined 
only on the integers. 


Solid straight line is vs. N 
representing the squares. The other 
curves are various analytic fits of the 
open circles to represent J The con- 
stant of proportionality has been chosen 
to be unity. 
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Table I: Selected values of -b r7 . 

NJ 



) 10(4) 

20(8) 

30(13) 

40(17) 

50(21) 

1 

262 

1123 

2584 

4646 

7307 

2 

187 

971 

2355 

4339 

6922 

3 

116 

822 

2128 

4035 

6541 

4 

47.4 

6 76 

1905 

3733 

6162 

5 


533 

1684 

3435 

5786 

6 


393 

1466 

3140 

5413 

7 


256 

1251 

2847 

5043 

8 


122 

1039 

2558 

4676 

9 



831 

2271 

4312 

10 



625 

1988 

3951 

11 



423 

1708 

3593 

12 



225 

1431 

3238 

13 



29-3 

1157 

2886 

14 




886-3 

2538 

15 




619.0 

2192 

16 




355-2 

1850 

17 




9.472 

1511 

18 





H75 

19 





842.7 

20 





513.7 

21 





188.2 


Table II: A comparison of the eigenvectors for N = 50 and 100 and for 
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Table III: Comparison of Energies (in ryd) . 


N { 

J D 

a 

e d 

j V v = R ) 

J w 

2 

1 

0.2 

-. 01579 a 

-.00407 

0 

3 

1 

0.2 

-.00992 

-.01272 

1 

4 

2 

0.19 

-.00409 

-.00996 

1 



0.09 

-.00012 

> 0 


5 

2 

0.17 

-.0025 

-.0074 






1 



0.11 

-.0005 

> 0 


6 

3 



-.0056 

_ p.— 





-.00030 






> 0 


...... . 




— ... ..... ... ... , ... 

...... 

7 

3 



-.0043 

s2 





-.00017 






> 0 


8 

3 



-.00345 

<s 2 





-.0014 






> 0 



-- — 


•—< * ■ 

— 



9 

4 



-.00291 



-.0014 

> 0 
> 0 


a. With symmetrization this value reduces to ^ - -0.0375 at a = 0.25 . 

b. This number is inferred by extrapolation from the values of . 
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Table IV: The effect of varying v in Y 

W 



N 

V 

lowest root 

second root 

8 

8 

-.003^5 

-.00141 


9 

-.00167 ! 

> 0 


10 

- . 00046 

> 0 


11 

+.00015 

> 0 

9 

9 

-.00291 

-.00141 


10 

-.00292 

-.00014 


11 

- . OOO98 

> 0 
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